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By using a model Hamiltonian with d-wave superconductivity and competing antiferromagnetic 
(AF) orders, the local staggered magnetization distribution due to nonmagnetic impurities in cuprate 
superconductors is investigated. From this, the net magnetic moment induced by a single or double 
impurities can be obtained. We show that the net moment induced by a single impurity corresponds 
to a local spin with S z =0, or 1/2 depending on the strength of the AF interaction and the impurity 
scattering. When two impurities are placed at the nearest neighboring sites, the net moment is 
always zero. For two unitary impurities at the next nearest neighboring sites, and at sites separated 
by a Cu-ion site, the induced net moment has S z =0, or 1/2, or 1. The consequence of these results 
on experiments will be discussed. 

PACS numbers: 75.20.Hr, 74.25.Jb, 72.10.Fk, 71.55.-i 



The nonmagnetic impurity effect in high temperature 
superconductors (HTS) has attracted significant interest 
both experimentally and theoretically for many years. 
The induction of local magnetic moment would be ex- 
pected due to the competition between spin magnetism 
and superconductivity in these systems. Nuclear mag- 
netic resonance (NMR) measurements in YBa2Cus07_ x 
have indicated that nonmagnetic Zn/Li impurities, en- 
hance antiferromagnetic correlation and staggered mag- 
netic moment is induced on the Cu ions in the vicinity of 
the impurity sites Q, 0, H 0, H E • 

Low-temperature 

scanning tunneling microscopy (STM) experiments Q 
have directly observed a sharp near zero bias resonance 
peak around the Zn impurity atoms on the surface of su- 
perconducting Bi2SroCaCu oO« 4-rr ( BSCCO). Both local 
potential scattering JlJI^jIl E3 13 

and Kondo im- 
purity scattering [13 . in <i-wave superconductors 
has been theoretically investigated. Recently Wang and 
Lee ^3 studied the formation of the local moment near a 
nonmagnetic impurity in HTS using a renormalized mean 
field theory of the t — J model. They tried to reconcile 
the formation of local moment around a unitary nonmag- 
netic impurity and strong zero bias resonance in the local 
density of states (LDOS) spectrum. 



We show that the net moment induced around the im- 
purity can be attributed to a local spin with S z =0, or 
1/2 depending on the impurity strength and the value 
of U, and we shall also present a phase diagram for the 
net moment formation. With multiple impurities, we ex- 
pect that quantum interference effect should exhibit it- 
self. When two impurities are placed at the n.n. sites, our 
numerical result indicates that the net induced moment 
is always zero, regardless of the value of U, the impurity 
strength, and doping. For two unitary impurities at the 
next n.n. sites and at sites separated by a Cu-ion site, 
we demonstrate that the net moment induced by them 
can be represented by a local spin with S z =0, or 1/2, 
or 1 depending on the strength of the AF interaction or 
doping. 

We begin with a phenomenological model Hamiltonian 
in a two-dimensional plane, in which both the DSC and 
the competing AF or the spin density wave (SDW) order 
are taken into account: 

H = - ^ tjj4 a c jo . + y^(Un id - + eS imti - ^)c\ a c i(T 

+ J2^ c h c ]i + h - c -^ 



On the other hand, quantum interference effect among 
multiple impurities in HTS draws considerable attention 
in recent years pj| Oil l2fl Liil UA, liij • The energy de- 
pendent of the spatial distribution of LDOS spectrum 
changes remarkably by varying the distance and orien- 
tations among the impurities. To our knowledge, so far 
there exists no study considering the local moment for- 
mation due to quantum interference effect among multi- 
ple impurities in HTS. In order to study this and related 
problems, we apply an effective model Hamiltonian de- 
fined on a square lattice with a nearest neighboring (n.n.) 
interaction Vdsc to simulate the d-wave superconductiv- 
ity (DSC) and an onsite Coulomb interaction U to repre- 
sent the competing antiferromagnetic (AF) order. In this 
paper, the local magnetic moment distribution induced 
around a single impurity will be numerically studied first. 



where e is the on-site impurity strength, i m is the im- 
purity site, jjL is the chemical potential. The hopping 
term includes the n.n. hopping t and the next n.n. hop- 
ping t r . The staggered magnetization and DSC order in 
cuprates are defined as Mf — (— l) l {c[^c^ — cj^c^) and 
\j = Vnsc(ci]Cji — C£|Cjf)/2. The mean-field Hamil- 
tonian (1) can be diagonalized by solving the resulting 
Bogoliubov-de Gennes equations self-consistently 
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where the single particle Hamiltonian Hij G = —Uj + 
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v ii u l^ ^ an k \ 2ksT J ' w ^ /(^) as tne Fermi distribu- 
tion function. The DSC order parameter at the zth site 
is Af = (Ai+ exyi + Ai_ eajJ i - A M+6y - A M _ e J/4 where 
denotes the unit vector along direction. Var- 

ious doping concentrations can be tuned by varying the 
chemical potential. In the present calculation, we set 
the lattice constant a and hopping integral t as units, 
t' = —0.2 and Vdsc = 1-0- Due to the localized na- 
ture of the impurity states, the order parameters around 
impurity are insensitive to the boundary conditions for 
large system sizes. The linear dimension of the unit cell 
of the vortex lattice is chosen as N x x N y — 32 x 32. The 
averaged electron density is fixed at n = 0.85. The calcu- 
lation is performed in very low temperature regime. The 
supercell techniques are employed to calculate the LDOS. 
The number of the unit cells is M x x M y = 20 x 20. 

Our numerical results for a single impurity show that 
the local AF order may be absent around the impurity 
site for small e and is present when impurity strength e 
becomes larger. In Fig. 1, we plot three typical spatial 
distributions of the staggered magnetization and their 
corresponding LDOS spectra around the impurity site. 
The net local moment is defined as S z = The 
impurity is situated at (16,16). Panels (a) and (b) cor- 
responds to U = 2.0 and e = 3. It is clear that no AF or- 
der induction is shown for such a weak impurity. In panel 
(b), the resonance energy at the impurity site (black line) 
is less than zero while the resonance energy at the n.n. 
impurity site (red line) is greater than zero. The green 
line is for the site far from the impurity site. With the 
increasing of the impurity strength, the resonance peak 
height at the impurity site becomes weaker and its peak 
shifts to zero energy while the resonance peak height at 
the n.n. impurity site turns stronger and its peak also 
moves to zero energy. Panels (c) and (d) are for a strong 
or unitary impurity e = 100 and Z7=2.0. We find the in- 
duced staggered magnetization reaches a maximum value 
of 0.08 at the n.n. site of the impurity in Fig. 1(c), and 
the local static AF fluctuation extends over several lattice 
sites from the impurity but there is no net induced mo- 
ment (or S z =0). The remarkable enhancement of zero 
bias resonance peak at the n.n. sites is shown in panel 
(d), reflecting the characteristics of a unitary impurity. 
The presence of weak local AF order results also a rather 
weak splitting of resonance peak. In Figs. 1(e) and 1(f) 
we show the our results for U = 2.35 and e = 100. A pro- 
nounced two-dimensional SDW order is clearly induced 
around this unitary impurity in Fig. 1(e). The staggered 
magnetization at the n.n. site reaches to 0.33 and the 
net induced moment becomes S z = 1/2. The zero bias 
resonance peak of the LDOS at a smaller U with S z = 
shown in Fig. 1(d) at the n.n. site is substantially sup- 
pressed and become two weak peaks in the larger U case 
with S z = 1/2 (see Fig. 1(e)). 

To examine the local moment formation, we present 
the phase diagram of e versus U in Fig. 2. It is obvious 
that the induced net moment (S z = 1/2) should show up 
for larger impurity strength e and stronger AF interaction 
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FIG. 1: Spatial variations of the staggered magnetization Mf 
[(a), (c) and (e)], and LDOS spectra [(b), (d) and (f)]. The 
impurity site is at (16,16). The black line in panels (b), (d) 
and (f) is for (16,16), red line for (17,16), blue line for (17,17) 
and green line for (16,1). The upper panels [(a), (b)], the 
central panels [(c), (d)] and the lower panels [(e), (f)] are for 
(U = 2.0, e = 3), (U = 2.0, e = 100) and (U = 2.35, e = 100), 
respectively. 

[/, while S z = tends to exist for smaller e or weaker U. 
In fact, there exist three phase regimes depending on the 
magnitude of U value or doping (not shown here). For 
small J7, no AF order (nonmagnetic phase) is induced 
around the impurity (see Fig. 1(a)). For intermediate U, 
weak local AF order (AF fluctuation phase) appears (see 
Fig. 1(c)). In both cases the net induced magnetic mo- 
ment corresponds to a local spin with S z = 0. For larger 
[/, a pronounced SDW order is induced (see Fig. 1(e)) 
and the net moment becomes S z — 1/2 (local moment 
phase). In the nonmagnetic phase regime and for a uni- 
tary impurity, there is a sharp zero bias resonance peak 
in the LDOS without any splitting 9]. The splitting is 
gradually showing up in the AF fluctuation phase regime. 
The LDOS could be substantially suppressed at zero bias 
and split into two weak peaks in the third regime with 
S z =l/2. It is important to notice that the AF fluctua- 
tion phase regime discussed here is absent in Ref. [17]. 
The staggered magnetization distribution around a uni- 
tary impurity like Zn obtained in this paper is in agree- 
ment with the results derived from the NMR experiments 
B II H II H H 13 The critical value of the U for induc- 
ing S z = 1/2 is doping dependent. The underdoped case 
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FIG. 2: Phase diagram of e versus interaction strength U for 
various phases near the impurity. 



corresponds to smaller U c . The S z = 1/2 moment due 
to a unitary impurity is much easier to be induced in 
an underdoped sample than in optimally and overdoped 
samples, because the amplitude of the induced AF or- 
der becomes increasingly stronger as the hole doping is 
decreased. Since the existence of imhomogenity in the 
HTS sample has been experimentally confirmed [2 ? 
as a result, the number of induced S z = 1/2 moments 
would be smaller than the number of unitary impurities 
like Zn. We also predict that the strong zero-bias-peak 
observed by Pan et al. |8| at the Zn impurity site should 
be associated with S z = moment in the overdoped and 
possibly optimally doped regions. In the underdoped re- 
gion, the induced moment would become S z = 1/2 where 
the LDOS spectrum should be much suppressed by the 
induced SDW at zero bias. 

We next study the quantum interference effect on the 
local moment formation due to two unitary impurities. 
When they are placed at the n.n. sites (see Fig. 3(a)), our 
numerical results for the distribution of the induced mag- 
netization around the impurities are shown in Figs. 3(b) 
and 3(c) with respectively 17=2.2 and 2.4. For 17=2.0 
and 2.35 (not shown here), the induced staggered magne- 
tizations are also uniformly zero and exactly identical to 
that in Fig. 3(b). Although it appears that the staggered 
magnetizations due to the two n.n. impurities have exact 
cancellation, their influence on the LDOS is still apparent 
and the result will not be present here. With [7=2.4 and 
no impurities, one can numerically demonstrate that the 
staggered magnetization has a stripe like structure [2(| 
with periodicity 8a which coexists with the DSC. The 
presence of the impurities could pin the stripes but does 
not modify the overall stripe-like structure except the 
magnetizations at or very close to the impurity sites are 
altered (see Fig. 1(c)). In all these cases, the net induced 
moment has S z = 0. This result is very robust and in- 
dependent of the value of £7, impurity strength e, and 
doping. For two unitary impurities placed at the next 
n.n. sites (see Fig. 3(d)), The induced spatial profiles 
of the staggered magnetization for U=2.0 and 2.4 are re- 
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FIG. 3: Impurities configuration [(a) and (d)] and staggered 
magnetization Mf [(b), (c) and (e), (f)]. The left panels [(a), 
(b), and (c)] and the right panels [(d), (e), and (f)] are for two 
n.n. impurities and two next n.n. impurities, respectively. 
Panels [(b), (c)] and [(e), (f)] correspond to U = 2.0 and 
U = 2.4, respectively. 



spectively shown in Figs. 3(e) and 3(f). The net moment 
associated with Fig. 3(e) yields a local spin of S z = 1/2, 
while that associated with Fig. 3(f) has S z = 1. For 
U = 2.35, the induced SDW no longer has the stripe like 
structure and we still obtain S z = 1. The induced net 
moment has been shown to have S z = when U is less 
than 1.9. 

Finally, we place two unitary impurities at sites sep- 
arated by a Cu ion (see Fig. 4(a)). The distributions 
of the staggered magnetizations are shown in Figs. 4(b), 
4(c) and 4(d) for {7=1.9, 2.0 and 2.4, and the induced 
net moments have respectively S z =0, 1/2 and 1. In Fig. 
4(b), no induced AF order is present for small U . The 
net moment has S z = 0. In Fig. 4(c), a remarkable 
enhancement of local AF order is shown at the central 
Cu site and the DSC order at this site is suppressed to 
almost zero. This is a kind of constructive interference 
effect. Its net moment has S z = 1/2. In Fig 4(d), one 
can clearly observe the pinning of SDW stripes by impu- 
rities. The local magnetization at the central Cu site is 
also enhanced. The net moment has S z = 1. As shown 
in Fig. 4(e), if we choose U = 2.2 and a different band 
parameter t r = —0.3, the induced local moment associ- 
ated with two individual impurities would have opposite 
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FIG. 4: Unitary impurities configuration (a) and spatial 
profiles of staggered magnetization Mf [(b), (c) ,(d) and 
(e)] in a 32 x 32 lattice. Their corresponding parameters 
are (b) (U = 1.9), (c) (U = 2.0), (d) (U = 2.4) and (e) 
(U = 2.2,t f = -0.3). 

polarity and yield a net S z = 0. The destructive inter- 
ference effect is shown at the center site where local AF 
order is equal to zero. This type of opposite polarities 
occurs when two impurities separated further apart even 
with t' = —0.2 as used in the present paper. The net 
local moment induced by two impurities spaced by one 
Cu ion could result S z = 0, or 1/2 or 1. In other words, 



quantum phase transitions may occur among different lo- 
cal moment phases by varying U values or doping. The 
detailed study of various configurations of the impurities 
on the local moment formation and LDOS spectra will 
be presented as a future work. 

In summary, we have investigated the induction of the 
local moment by a single and double impurities in HTS 
based on a phenomenological model with DSC and com- 
peting AF orders. By tuning the impurity potential and 
the value of U, a transition between various net magnetic 
moment states may appear. The LDOS has been calcu- 
lated for a single impurity. We show that the zero bias 
resonant peak obtained next to the unitary impurity site 
is always associated with weak U and S z = 0. When 
S z = 1/2, the LDOS at zero bias is suppressed by the 
gap of the locally induced SDW. This is consistent with 
the recent STM experiments 25], where the zero bias 
resonant peaks due to Zn impurities are observed only 
in the hole rich region, not in the hole poor region in 
BSCCO. In addition, the quantum interference effect by 
two nonmagnetic impurities has also been studied. Our 
calculation predicts the absence of net magnetic moment 
around two n.n. impurities, regardless of the values of 
the impurity strength, doping, and U. This result indi- 
cates that the number of induced S z = 1/2 moments is 
always smaller than the number Zn impurities even in an 
underdoped sample where the AF strength is apprecia- 
ble. We demonstrate that net magnetic moment around 
two unitary impurities placed at the next n.n. sites and 
sites separated by a Cu ion could be either 0, 1/2 or 
1 depending on the value of U. It is our hope that the 
present investigation on the local moment formation may 
provide useful information for future experimental test. 
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